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ABSTRACT 



Using Particle-In-Cell simulations i.e. in the kinetic 
plasma description the discovery of a new mechanism 
of parallel electric field generation was recently reported. 
Here we show that the electric field generation parallel to 
the uniform unperturbed magnetic field can be obtained 
in a much simpler framework using the ideal magnetohy- 
drodynamics (MHD) description. In ideal MHD the elec- 
tric field parallel to the uniform unperturbed magnetic 
field appears due to fast magnetosonic waves which are 
generated by the interaction of weakly non-linear Alfven 
waves with the transverse density inhomogeneity. In the 
context of the coronal heating problem a new two stage 
mechanism of plasma heating is presented by putting em- 
phasis, first, on the generation of parallel electric fields 
within an ideal MHD description directly, rather than fo- 
cusing on the enhanced dissipation mechanisms of the 
Alfven waves and, second, dissipation of these parallel 
electric fields via kinetic effects. It is shown that for a sin- 
gle Alfven wave harmonic with frequency 1/ = 7 Hz, and 
longitudinal wavelength Xa = 0.63 Mm for a putative 
Alfven speed of 4328 km s ^ ^ , the generated parallel elec- 
tric field could account for 10% of the necessary coronal 
heating requirement. We conjecture that wide spectrum 
(10^^ — 10'' Hz) Alfven waves, based on the observa- 
tionally constrained spectrum, could provide the neces- 
sary coronal heating requirement. By comparing MHD 
versus kinetic results we also show that there is a clear 
indication of the anomalous resistivity which is 100s of 
times greater than the classical Braginskii value. 

Key words: Sun: oscillations - Sun: Corona - (Sun:) 
solar wind. 



1. INTRODUCTION AND MOTIVATION 



The coronal heating problem, the puzzle of what main- 
tains the solar corona 200 times hotter than the photo- 
sphere, is one of the main outstanding questions in so- 
lar physics. A significant amount of work has been done 
in the context of heating of open magnetic structures in 
the solar corona (e.g. phase mixing, one of the possi- 
ble mechanisms of the heating). Historically all phase 



mixing studies have been performed in the Magnetohy- 
drodynamic (MHD) approximation, however, since the 
transverse scales in the Alfven wave collapse progres- 
sively to zero, the MHD approximation is inevitably vi- 
olated. Thus, Tsiklauri et al. (2005a b) studied the phase 
mixing effect in the kinetic regime, using Particle-In-Cell 
simulations, i.e. beyond a MHD approximation, where 
a new mechanism for the acceleration of electrons due 
to the generation of a parallel electric field in the solar 
coronal context was discovered. This mechanism has 
important implications for various space and laboratory 
plasmas, e.g. the coronal heating problem and acceler- 
ation of the solar wind. It turns out that in the magneto- 
spheric context, a similar parallel electric field generation 
mechanism in the transversel y inhomogeneous plasma s 
has been previously reported JGenot et all 120041 1 19991) . 
See also Mottez et al. ( 200d) and references therein. 

After the comment paper by 'Mott ez et alJ ([2006') we 
came to the realisation that the electron acceleration 
seen in both series of works jTsiklauri et all l2005albt 
iGenot et'aill2004ll999 h is a non-resonant wave-particle 
interaction effect. In works by Tsiklauri et al. (2005a b) 
the electron thermal speed was vth,e — 0.1c while the 
Alfven speed in the strongest density gradient regions 
was va = 0.16c; this unfortunate coincidence led us to 
the conclusion that the electron acceleration by parallel 
electric fields was affected by the Landau resonance with 
the phase-mixed Alfven wave. In works bv iGenot et aTl 
( I2OO4. 1999.) the electron thermal speed was Vth,e = 0.1c 
while the Alfven speed was va — 0.4c because they 
considered a more strongly magnetised plasma applica- 
ble to Earth magnetospheric conditions. However, the 
interaction of the Alfven wave with a transverse den- 
sity plasma inhomogeneity when the Landau resonance 
condition uj = k^^CA{x) is met can be quite i mpor- 
tant for the electron acceleration (Chaston et al., ■20001 
[Hasegawa & Chen, 1976). Chaston et al. (2000) assert 
that the electron acceleration observed in density cavi- 
ties in aurorae can be explained by the Landau resonance 
of the cold ionospheric e lectrons with the Alfven wave. 
iHasegawa & Cher ? (1976) also established that at the res- 
onance Alfven wave fully converts into the kinetic Alfven 
wave with the perpendicular wavelength comparable to 
the ion gyro-radius. We can then conjecture that because 
of kinetic Alfven wave front stretching (due to phase mix- 
ing, i.e. due to the differences in local Alfven speed). 



this perpendicular component gradually realigns with the 
ambient magnetic field and hence creates the time vary- 
ing parallel electric field component. This points to the 
importance of the Landau resonance for electron accel- 
eration when the reson ance condition is m et. But as 
witnessed from works of Genot et al. (2004, 1999), even 
when the resonance condition is not met electron acceler- 
ation is stiU possible. 

There were three main stages that lead to the formulation 
of the present model. 

(i) The realisation of the parallel electric field generation 
(and particle acceleration) being a non-resonant wave- 
particle interaction effect lead us to the question: could 
such parallel electric fields be generated in a MHD ap- 
proximation? 

(ii) Next we reahsed that if one considers non-linear gen- 
eration of the fast magnetosonic waves in the transversely 
inhomogeneous plasma, then E = ~{V x B)/c contains 
a non-zero component parallel to the ambient magnetic 
field = -{V^By - VyB,)/c. 

(iii) From p revious studies (iBotha et aLl l200rt 
iTsiMauri et all l2001l) we knew that the fast magne- 
tosonic waves (Vx and B^) did not grow to a substantial 
fraction of the Alfven wave amplitude. However after re- 
producing the old parameter regime (k=l, i.e a frequency 
of 0.7 Hz), the case of k=10, i.e a frequency of 7 Hz was 
considered, which showed that fast magnetosonic waves 
and in turn a parallel electric field were more efficiently 
generated. 



2. MODEL, RATIONALE, AND MAIN RESULTS 

Unlike pre vious studies JTsiklauri et all l2005iibl: 
iGenot et all |2004, 199§, here we use an ideal MHD 
description of the problem. We solve numerically ideal, 
2.5D, MHD equations in Cartesian coordinates, with 
a plasma beta of 0.0001 starting from the following 
equilibrium configuration: A uniform magnetic field Bq 
in the z— direction penetrates plasma with the density 
inhomogeneity across the a;— direction, which varies 
according to 

p{x) = po [1 + 2 (tanh(a; + 10) + tanh(-.T + 10))] . 

(1) 

This means that the plasma density increases from some 
reference background value of pQ, which in our case was 
fixed at po = 2 X lO^prrip g cm^'^ (with a molecular 
weight of ^ — 1.27 corresponding to the solar coronal 
conditions ^H:^He=10:l and rrip being the proton mass), 
to 5po- Such a density profile across the magnetic field 
has steep gradients with a half-width of 3 Mm around 
X ~ ±10 Mm and is essentially flat elsewhere. Such 
a structure mimics e.g. the footpoint of a large curva- 
ture radius solar coronal loop or a polar region plume 
with the ratio of the density inhomogeneity scale and the 
loop/plume radius of 0.3, which is the median value of 
the observed range 0. 15 - 0.5. 



The initial conditions for the numerical simulation are 
By — v4cos(fcz) and Vy = —CA{x)By at t = 0, which 
means that a purely Alfvenic, linearly polarised, plane 
wave is launched travelling in the direction of positive 
zs. The rest of the physical quantities, Vx and Bx (which 
would be components of fast magnetosonic waves if the 
medium were totally homogeneous) and Vz and B^ (the 
analogs of slow magnetosonic waves) are initially set to 
zero. The plasma temperature is varied as the inverse 
of Eq.(l) so that the total pressure always remains con- 
stant. Boundary conditions used in our simulations are 
periodic along the z- and the zero gradient along the x- 
coordinates. We fixed the amplitude of Alfven wave A 
at 0.05 throughout. This choice makes the Alfven wave 
weakly non-linear. 

As a self-consistency test, we considered situation when 
the wavenumber of the initial Alfven wave is A; = 1. 
In dimensional units this corresponds to an Alfven wave 
with frequency (i/ — 0.7 Hz), i.e. longitudinal wave- 
numbers A4 = 6 .3 Mm. We corroborated the previous 
results of lBodia et al. (2000) and Tsiklauri et al. (.2001) . 

We now present results for the case of large wave- 
numbers, fc = 10, which in dimensional units correspond 
to an Alfven wave with v = 1 Hz, and = 0.63 Mm. 
This is a regime not investigated before. In Fig. (1) we 
show shaded surface plots of both the fast magnetosonic 
waves (14) and parallel electric field E^_. We gather from 
this gra ph th at similarly t o the r esults of Tsiklauri et alj 
(12005b') and iGenot et al. (2004), the generated pai-allel 
electric field is quite spiky, but more importantly large 
wave-numbers i.e. short wavelength now are able to sig- 
nificantly increase the amplitudes of both the fast mag- 
netosonic waves (14) and parallel electric field E^- This 
amplitude growth is beyond a simple scaling. The am- 
plitude growth is presented quantitatively in Fig. (2) (by 
doubling spatial resolution we perform satisfactory con- 
vergence test). The amplitude of t4 now attains values 
of 0.01 unlike for moderate fcs. Thus, large wavenumbers 
(i.e. stronger spatial gradients) seem to yield larger val- 
ues for the level of saturation of the I4 amplitude. In the 
considered case, E^ now attains values of 0.001. 

The coronal energy losses that need to be compensated 
by some additional energy input, to keep the solar corona 
at the observed temperatures, are (in units of erg cm^^ 
s^^): 3 X 10^ for the quiet Sun, 8 x 10^ for a coronal 
hole and 10^ for an active region. One can estimate the 
heating flux per unit area (i.e. in erg cm^^ s^^): 

^.^i..A.^5xl03(^)^(^), (2) 

where Eh ~ 10^^ erg cm^"^ s^^. This yields an estimate 
of Fh « 2 X 10^ erg cm~-^ s^^ in an active region with a 
typical loop base electron number density of — 2x 10^ 
cm-3 and T 1 MK. 

The energy density associated with the parallel electric 
field Ez is 

Ee^^, [ergcm-3] (3) 







Figure 1. Snapshots ofVx and at t — 2 and 20 for the case of k = 10, v = 7 Hz, Xa = 0.63 Mm. 



where e is the dielectric permitivity of plasma. The latter 
can be deduced from 



Anpc 
~B2- 



(4) 



For the coronal conditions (p = 2 x 10^/imp g cm ^, 
p = 1.27, B = 100 Gauss) e « 4.8048 x 10^. In Fig. (2) 
we saw that electric field amplitude attains a value of 
w 0.001. In order to convert this to dimensional units 
we use = 4328 km s^^ and B = 100 G and E^, = 
-{V^By -VyBa:)/cto obtain E:, « {c°^B/c) x 0.001 = 
0.0014 statvolt cm^^ (in Gaussian units). Therefore the 
energy density associated with the parallel electric field 
Ez (From Eq.(3)) is 



EE = ex 0.0014V(87r) = 3.7471 x 10" 



(5) 



[erg cm 



-31 



In order to get the heating flux per unit area for a sin- 
gle harmonic with frequency 7 Hz, we multiply the latter 
expression by the Alfven speed of 4328 km s^^ to obtain 



Eec\ 



1.62x10^ [ergcm-^s"^] 



(6) 



which is « 10 % of the coronal heating requirement es- 
timate for the same parameters made above using Eq.(2). 
Note that the latter estimate is for a single harmonic with 
frequency 7 Hz. 

A crucial next step that is needed to understand how the 
generated electric fields parallel to the uniform unper- 
turbed magnetic field dissipate must invoke kinetic effects. 
In our two stage model, in the first stage bulk MHD mo- 
tions (waves) generate the parallel electric fields, which 
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Figure 2. Time evolution of the amplitudes of Vx = 
and Ez = Solid lines with stars represent solu- 

tions using the Lare2d code with 4000 x 4000 resolution, 
while dash-dotted lines with open symbols are the same 
but with 2000 x 2000 resolution. Here A = 0.05, k = 10, 
v = 1Hz,\a = 0.63 Mm. 



cannot accelerate particles if we desc ribe plasmas in the 
ideal MHD limit. Genot et al. (2004') andl Tsiklauri et alJ 
(2005b) showed that when the identical system is mod- 
elled in the kinetic regime particles are accelerated with 
such parallel fields and Alfven wave energy is converted 
into heat on a time scale of a few Alfven periods. 

Alfven waves as observed in situ in the solar wind al- 
ways appear to be propagating away from the Sun and 
it is therefore natural to assume a solar origin for these 
fluctuations. However, the precise origin in the so- 
lar atmosphere of the hypothetical source spectrum for 



Alfven waves (turbulence) is unknown, given the impos- 
sibility of remote magnetic field observations above the 
chromosphere-corona transition region. Studies of ion 
cyclotron resonance heating of the solar corona and high 
speed winds exist which provide important spectroscopic 
constraints on the Alfven wave spectrum. Although the 
spectrum can and is observed at distances of 0.3 AU, it 
can be projected back to the base of corona using empiri- 
cal constraints. Therefore, we conjecture that wide spec- 
trum (10^^ — 10'' Hz) Alfven waves, based on the obser- 
vationally constrained spectrum, could provide the nec- 
essary coronal heating requirement. The exact amount of 
energy that could be deposited by such waves through our 
mechanism of parallel electric field generation can only 
be calculated once a more complete parametric study is 
done. Thus, the "theoretical spectrum" of the energy 
stored in parallel electric fields versus frequency needs 
to be obtained. At present we only have two points, 0.7 
Hz and 7 Hz, in our "theoretical spectrum". Preliminary 
results will be presented elsewhere cTsiklauri ■■2006a-b.) . 



dVA/dy we put w VA{y*)/Sy where Va(?/*) is the 
Alfven speed at the strongest gradient point and l/(5y is 
the strength of the gradient. The latter can be approxi- 
mated as (Va(oo) - VA{0))St/{50S) ■ 1/(50(5), with 50S 
being the scale of the transverse density gradient. Putting 
u! — O.Swci and 6t = 54.69/wcj we obtain ri = 7 x 
10_fVl(y*)A/(VA(cx))-VA(0))2. From Tsiklauri et alJ 
iHE) VAiy*) = 0.16c and {Va{^) - ^^(0)) = 
0.25c - 0.125c = 0.125c, hence r/ = 1.8 x lO'^cA. 
Also, A ~ vth^e/'-^ce- The thermal speed of electrons (at 
infinity) Vth,e = \/kT/mQ, which for IMK coronal tem- 
perature is 3.9 X 10^ cm s^^. The electron cyclotron fre- 
quency ijjce = 1-76 X 10^i3[G'], which for 10 G magnetic 
field is 1.8 X 10* s^^. Thus, A = 2 cm and we finally ob- 
tain 77 = 1.1 X 10^ cm-^ s^^, or in SI units 77 = 1.1 x 10^ 
s^^. This is about 100 times greater than the clas- 
sical Braginskii value for the resistivity. Therefore, the 
obtained value is a clear indication of the anomalous re- 
si stivity. Inter estingly other s quote simi lar values for the 
r? jPetkaki et al.i.2003; .Wattltlill2(X)l . 



3. MHD VERSUS KINETIC EFFECTS AND 
ANOMALOUS RESISTIVITY 



There are two main candidates for the solution of the 
coronal heating problem: so-called DC and AC models. 
DC or magnetic reconnection based models need to in- 
voke anomalous resistivity (somewhat ad hoc concept, 
but there is a lot of indirect evidence for it). At present the 
details of reconnection in 3D are not understood in full. 
AC models are often not quantitative, or not enough heat- 
ing can be provided, unless some enhanced dissipation 
mechanisms (e.g. phase-mixing, resonant absorption) are 
invoked. This makes a perfect playground for studying 
interplay between MHD and kinetic theories of the solar 
corona. The reason is two-fold: 1) The issue of anoma- 
lous resistivity on which DC models rely can only be set- 
tled by studying kinetic effects (micro-physics); 2) AC 
models which use MHD eventually break down, as of- 
ten system naturally evolves towards progressively small 
scales (e.g. in phase-mixing); Therefore, kinetic effects 
become important. 



iTsiklauri et al.l ll2005albh . amongst other findings, es- 
tablished that the Alfven wave amplitude decay law in 
the inhomogeneous regions, in the kinetic regime is oc 

exp - (a;A/1250A)-' (where X is the coordinate along 

uniform magnetic field - mind the change of the geom- 
etry!); which is the same as in the MHD approxima- 
tion discovered by Heyvaerts and Priest (1983) (HP83 



thereafter): oc cxp 



6V? 



Question that 



begs to be asked: what if we calculate resistivity from 
the latter MHD formula using our kinetic (PIC) em- 
pirical dissipation length of 1250A? (F. Malara, pri- 
vate communication). By equalising expressions un- 
der the both exponents and noting that x from HPS 3 
is equal to xA from our kinetic formula, we obtain 
■q = 6y|/[w2v^2(l250A)3]. In order to estimate = 
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